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1. Introduction 

Recent observations of high redshift supernovae suggest that the universe is currently 
undergoing a phase of acceleration ||l|, Q. This is usually explained by the presence 
of 'dark energy', some as yet unknown, negative pressure fluid. In the standard 
model of cosmology, the ACDM model, this dark energy is assumed to be vacuum 
energy in the form of a small, positive cosmological constant. Measurements of the 
cosmic microwave background (CMB) anisotropics by WMAP and large scale 
structure surveys [Q suggest that this fluid must make up ~ 70% of the content of 
the universe, the remainder being made up of matter, both baryonic matter (~ 4%) 
and dark matter. Despite the fact that the ACDM model is currently the best fit 
to the cosmological data, there are theoretical issues with the cosmological constant. 
For a naive estimate of A, there is a difference of around 120 orders of magnitude 
between the theoretical prediction and the value inferred from observations. As yet, 
no satisfactory explanation has been put forward to resolve this discrepancy and 
explain the observed smallness of A. 

Rather than trying to explain why the cosmological constant is so small, vari- 
ous alternative mechanisms for achieving late time acceleration have been explored. 
These can essentially be divided into models which modify the matter content of 
the universe, such as quintessence P|, or the more phenomenological Cardassian 
and Hobbit models 0, and those which modify the gravitational interaction, such 
as MOND [^, f{R) theories [^, and braneworld models. The braneworld scenario 
is a set-up in which we have extra dimensions in nature, which are hidden because 
we are confined to live on a slice - a brane in the higher dimensional spacetime. 
The most common models, ADD |TD| and Randall-Sundrum (RS) |TT| , lead to an 



- 1 - 



effective theory of braneworld gravity which is Einstein gravity at large scales, but 
with small scale Kaluza-Klein modifications. However, very early on in braneworld 
research, it was realised that braneworlds could also display large scale modifications 
of gravity [12, 13, 14, 15|. The DGP model in particular |1^ has received a great 



deal of attention as a possible viable cosmological alternative to ACDM. However, 
while being attractive from the phenomenological point of view, the DGP model has 
inconsistencies, such as ghosts JT^ [T^, pressure singularities and tunnelling instabil- 



ities ||T9|, (although see ||21| for counter-arguments). As such, it is surprising that 
so much focus has centered on DGP in comparison to other braneworld modified 
gravities. 

DGP, like many braneworld models, has a Z2 symmetry around the braneworld, 
but interestingly if one relaxes this symmetry, it is also possible to get IR modifica- 



tions of gravity such as the asymmetric models of Padilla, [^2|, 23|, and the hybrid 
asymmetric-DGP type "stealth" model of Charmousis, Gregory and Padilla (see 
also [^]). In each of these models, the cosmological constant and the Planck masses 
are different on each side of the brane. In the asymmetric model (which we focus on 
here) there is a strong hierarchy between the Planck masses and the adS curvature 
scales on each side of the brane. On one side of the brane, there is a large cosmo- 
logical constant and Planck mass and the interior of the bulk is retained, on the 
other side of the brane the cosmological constant and Planck mass are low, and the 
exterior of the bulk is kept. Keeping the bulk interior produces a localizing effect on 
the braneworld gravity, whereas keeping the exterior tends to produce an opening 
up effect and modifies gravity in the infrared [^]. With a judicious choice of scales, 
it is possible to have a regime in which gravity is effectively 4D, before opening up 
at very large scales. (At small scales of course, the KK modes cause gravity to be 
effectively 5D.) In ^ this model has been extensively tested from the point of 
view of particle physics, but it is less clear that it will pass cosmological tests and in 
particular, reduce to standard 4-dimensional General Relativity at early times. 

In this paper we explore the cosmology of the asymmetric (AC) model, focussing 



in particular on the type la supernova data and the expansion parameters from 
WMAP p, We first show how the AC model can be viewed as a one parameter 
extension of the DGP model over a wide range of scales. We then explore the effect of 
this additional parameter on the expansion history of the universe, making reference 
to the Supernova and WMAP data. Finally, we discuss how the additional parameter 
makes it more difficult to simultaneously fit the various observational constraints. We 
also comment on the inclusion of a bulk black hole. 



2. Asymmetric Braneworld Models 

We start by reviewing the asymmetric branes model, and deriving the cosmological 
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equations. The action can be written as |23] 



S — Shulk ~\~ Sfjrane (2-1) 

where the 'bulk' action contains the gravitational field dynamics, and is given by the 
Einstein Hilbert and Gibbons Hawking terms: 



2 brane 



x^f^K^''^. (2.2) 



Here, gab is the bulk metric with corresponding Ricci scalar R. The metric induced 
on the brane is 

lab = gab - naflb (2.3) 

where n"" is the unit normal to the brane in A4i pointing out of A^j. The extrinsic 
curvature of the brane in A^j is given by 

^i?=7:7.'V(,n,). (2.4) 

The brane action for the asymmetric branes model is given by 

j4. 



Sbrane — d x(— CT a/— 7 + C matter) (2-5) 

J brane 

where a is the brane tension and Cmatter describes the matter content on the brane. 

Note that with braneworld models, there is a subtlety in how we encode the grav- 
itational action. We can either view the brane as a genuine zero thickness object - a 
mathematical boundary between two different spacetimes (which just happen to be 
mirror images of each other in the usual Z2 braneworlds, such as Randall-Sundrum) 
- or as the zero thickness limit of some finite size object, an approximation to the 



domain wall of the early braneworld models PD|. The equivalence of these two de 



scriptions has been well established pO|, as well as an understanding of the next 



to leading order corrections (Although see for an interesting discussion of 
possible cosmological consequences of finite width.) These different physical perspec- 
tives, boundary vs. ^-function, translate into a different way of expressing the action, 
i.e. whether or not we use the Gibbons Hawking term. In the asymmetric case, since 
the Planck mass is different on each side of the brane, the boundary description 
is somewhat more natural, and makes it easier to obtain the correct equations of 
motion. 

The equations of motion for each bulk are given by the Einstein equations 

Rab - ^Rgab = -K9ab (2.6) 



while the brane equations of motion are found from the Israel conditions |^D[ to be 

2{M^Kab) - ^lab = \ {Tab ~ ^T^ab) , (2.7) 



3 



obtained by varying (|2.1|) with respect to the induced brane metric ^ab- The energy- 
momentum tensor for the additional matter on the brane is given by 

Tab = -^f^. (2.8) 

The background metric (jab is found by solving the equations of motion with 
Tab = 0, and may be written as 

ds'^ = gabdx°'dx'' = a^{y)rj^ydx^dx" + dy"^ (2.9) 

where = (x^, y) with the brane at y = 0, and a{y) is the warp factor which has 
the general form: 

a^{y) = e-'^^^^^y^ (2.10) 

where 9i = ±1, the subscript i = 1,2 refers to the two sides of the brane {i = 1 being 
y < 0), and Aj = —6kf defines the adS curvature scale on each side. 



The metric ( p. 91) and the equations of motion also impose a condition on the 
brane tension: 

2{Mfe,h) = ^ (2.11) 
o 

where {Z) = {Zi + Z2)/2 and AZ = Zi — Z2 for a quantity Zj differing across the 
brane. 



Three separate cases of this model were considered in |]23| for different 6i values: 
(i) the Randall-Sundrum (RS) case, 9i = 62 = I, (ii) the inverse Randall-Sundrum 
case, 02 = Oi = —1, and (iii) the mixed case, where 9i = —62 = 1. If Oi corresponds 
to the left-hand side of the brane [y < 0) and 62 corresponds to the right-hand side 
of the brane (?/ > 0), then the RS case has the warp factor decaying away from the 
brane on both sides, while the inverse RS case has the warp factor growing on both 
sides. In the mixed case, the warp factor decays away from the brane on the left. 



whilst growing on the right. As explained in detail in P3| , |26[| , whereas 4-dimensional 
Einstein gravity cannot be reproduced at any scale in the inverse RS case, it can 
be achieved in the RS and mixed cases, along with infra-red (IR) modifications. 
However, only the mixed case will approach a de Sitter state at late times, leading 



to exponential late-time acceleration without an effective cosmological constant [p^ ]. 
Therefore, only the mixed case where 9i = —62 = 1 is considered from now on. 

Turning to cosmological solutions, since we have Einstein gravity in the bulk, we 



know that the bulk is completely specified by the AdS-Schwarzschild metric [33 



where 



dv'^ 

-hi{r)dt^ + — - + r'^d^l (2.12) 



h,{r)=Pk^ + ^,-f^ . (2.13) 
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For simplicity, we will take the case where there is no black hole in either bulk, 
/ij = 0, and also consider only flat f2 = 1 universes, k = 0, hence 



hi{r) = r kj^ . 



(2.14) 



In order to construct the brane, we glue a solution in A4i to a solution in A^2, 
where the brane will form the common boundary. Then, in Aii, the boundary dAii 
is given by the section (^^(t), ai(r), x'^) of the bulk metric (|2.12|) , where r is the 
proper time of an observer comoving with the boundary, so that 



hi{ai) 



-1 



(2.15) 



where the differentiation is with respect to r. The outward pointing unit normal to 
dA4i is now given by 

na = ei{~a,{T),U{T),0) (2.16) 

where 6i = ±1 as before. For 6i = 1, Aii corresponds to < a < aj(r), while for 
6i = —1, A4i corresponds to aj(r) < a < cxo. The induced metric on dA4i is that of 
a FRW universe, 

ds"^ = -dr^ + aldyil. (2.17) 

Since the brane coincides with both boundaries, the metric on the brane is only well 
defined when ai(r) = 02 (r) = a(r) and the Hubble parameter is then defined as 
H = ^. If we now introduce a homogeneous and isotropic fluid on the brane, whose 
energy-momentum tensor is given by |Q 

(2.18) 



Tab = {p + P)Tan+V1a 



with energy density p, pressure p and r°, the velocity of a comoving observer (which 
in A4i is = {tiij), d(r), 0)), we can evaluate the spatial components of (pT7|). Doing 
this and using ( |2.15| ) to substitute for t, we find 
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(2.19) 



Substituting for a using ( p.ll| ), and h{a) using (p.l4|) , the modified Friedmann equa- 
tion for the mixed case is 



p = 6 



Mfym+kj - k,^ - Ml + ki - h 



(2.20) 



From ||2^, |26|, we know that there is a range of scales over which gravity is four 
dimensional, given by 



k^^ <^ r <^ rc 



Mf 
Miki 



(2.21) 
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which clearly requires Mi ^ M2. For this model to be phenomenologically viable, 
this range of scales must be appropriate. Since we are looking at Tc as representing 
the scale at which late time acceleration sets in, we expect the crossover scale to be 
of order the current horizon size, ~ Hq^. On the other hand, table-top tests of 
General Relativity have confirmed its validity down to sub-mm scales, which fixes 
our largest frequency scale, ki (the UV cut-off of the theory), so that ^ ~ 0.1mm. 
These constraints give us a large hierarchy of scales, and, as already noted, require 
a large hierarchy in the parameters. 

It is interesting to see these scales emerge from an analysis of the Friedmann 
equation (|2.20|) . Obviously ( p.20|) looks nothing like the standard Friedmann equa- 



tion, and so can only reduce to such in certain asymptotic limits. Since 

^— (k + ^ hTH<t:k 

\^IPT¥ ~ < 2fc 2.22 
[H for if > A; 

we see that we can only get the behaviour required if ki ^ ^2, and we take 
ii -C /ci. In this regime, the Friedmann equation can be written as 

Mf 



P 



3-^ii2 - 6Ml (^^H^ + kl - k^ (2.23) 



We therefore see the existence of an accelerating vacuum whenever 

Hl = ^k,^AkM-k^>^ (2.24) 



Ml \ Ml 

We can also read off the 4D Planck mass = l/SnG by comparing with the 
standard 4D Friedmann equation: 

p = SmliH' (2.25) 



as 

~ > (2.26) 



This agrees with the expression derived in and also with a direct computation 
of the propagator (see appendix 

We would like to compare ( |2.23| ) with the cosmological equations from the DGP 
model. The DGP model is characterized by an induced curvature term on the brane, 
and (in its original form) is Z2 symmetric around the brane, which is tensionless and 
embedded in 5D Minkowski space [|l4| : 

5dgp = M| j £xVGR(5) + Ml j d^xv^i?(4) (2.27) 
The brane cosmological equations from this action are given by |jT6| 



p = GM^H^ T 12MlH (2.28) 
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The choice of sign in the hnear Hubble term is due to the choice of which part of 
the bulk is kept. The minus sign, corresponding to the exterior being kept, is the 
self-accelerating hianch, which has late time cosmological acceleration. The crossover 
scale tdgp = Ml/2M^ corresponds to the scale at which gravity ceases to be 4D, 
and the extra dimension opens up. 

To compare the asymmetric and DGP models, note that if we take H ^ k2, 
then we may approximate the second bracket in ([2. 231), and obtain 



(2.29) 



which is of course ( |2.28|) after suitable substitution. Therefore, over a large range of 
scales, asymmetric cosmology can be viewed as a generalization of DGP cosmology. 
To parametrize this in a simple way for our analysis, we set 



a 



Ho Ml 



(2.30) 



where Hq is the current value of the Hubble parameter. ( p.23| ) then becomes: 



P 



^mliHl 



(2.31) 



Here, a = {2HQr ncp)'^ is essentially the same as the DGP crossover scale, and j3 is 
the new parameter coming from the asymmetric physics. It is precisely the effect of 
this new parameter which we seek to explore. 



3. Asymmetric Cosmology 

In order to explore the effect of the AC model, it is useful to rewrite the Friedmann 
equation in an Einstein form by solving ( p. 31 ) for E = H/Hq : 



E(z 



(zf = ni{l + 2)3(i+-0 + 2a{a -P) + 2a-^(a-/3)2 + a(l + ;2)3a+-0. (3.I) 

Here, an implicit sum over the various contributions to the energy density with 
equations of state pi = WiPi is understood. Note that the + root of the quadratic is 
required to get the correct — limit of the Israel equations. We can now readily 
compare the AC model with ACDM and DGP, which are implicitly contained in 
(13): a = and we include an for ACDM, and /3 = for DGP. Since the DGP 
model has been carefully analysed with cosmological expansion data (see e.g. pSi), 
here we focus qualitatively on the additional features the /3-term brings relative to 
DGP. 

The aim of gravitationally driven late time acceleration is to avoid using a cosmo- 
logical constant (I^a), therefore evaluating ( p.l|) at the current time gives a constraint 
between the model parameters a, (3 and the current matter density: 

nm=l-2a{^/TT^-(3) (3.2) 
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(ignoring the relatively insignificant radiation component). This means that once Qm 
is fixed, the asymmetric cosmology forms a one parameter family of solutions (note 
that DGP is entirely constrained by fixing 0,^ in a fiat universe). From ( |2.30| ), we 



see that both a and j3 are positive, and in addition self acceleration requires P < a 
from ( |2.24| ). Thus our additional parametric degree of freedom in the asymmetric 



model has a fairly limited range. 

The modified Friedman equation (|3.1|) shows clearly the effect of f3 over the 



range (0, a). As already noted, (3 = corresponds to the DGP model, with = 
in the usual notation of encoding the DGP crossover scale as an effective DGP Q 
contribution. The other limit, j3 = a, corresponds to an n = 1/2 Cardassian model 
0, or, alternatively, a dark energy fiuid with (constant) equation of state w = —1/2. 

As with DGP, relaxing the constraint of fiatness leads to a wider range of pa- 
rameter choice: 

(1 — — Qrn) ^ ^ (1 — ^A: — ^m) oN 

> a > -^=^= — (3.3) 



2^/rTj^ 2^1 — 

and the cosmological models now form a three parameter family, which can be la- 
belled using {^m, ^k, C(^} (or by trading one of the fl parameters for /?). In order to 
more readily compare with DGP results, we will use the former parametrization, and 
compare results in the {^Im, plane (recall = ilr^ i^i DGP) either for various 
fixed P values with Qk varying, or fixed flk values with /? varying. In spite of the 
enlarged parameter space, the asymmetric cosmology turns out to be under more 
cosmological tension than DGP. 

A nice way to encode this information is to consider the effective dark energy 
which is the difference between the square of the Hubble parameter and the matter 



content |]36[|: 



Qde{z) = 2a{a - p) + 2a ^ {a - (3^ + ^(1 + 2)3a+-«) , (3.4) 

and the effective dark energy pressure (again, the discrepancy between the Einstein 
pressure and the actual pressure): 



^de{z) = -a 



Using these, we can find an effective equation of state, 

(l+«;i)fii(l + z)3(i+»') 



(3.5) 



wde{z) = -1+ 



2[{a - (3)^/{a - (3^ + ^(1 + z)3(i+«'.) + (a - (3^ + ^(1 + z)3(i+'"»)] 

(3.6) 

This shows how the equation of state always has w > —1 for a > (3, and thus the 
model can never enter a phantom regime. We also see how for j3 > 0,w is raised from 
its DGP {(3 = 0) value (see figure |I|). Overall therefore, we expect that expansion 
data will favour a lower in both DGP and asymmetric models. 
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Figure 1: Variation of the effective dark energy equation of state taking = 0.25. Tlie 
values of a run from the minimum ahowed (/? = 0, the DGP value) in purple, to the 
maximum value, a = /?, shown in red. These clearly show how increasing /3 neutralizes 
acceleration in asymmetric cosmology. 



In order to see this explicitly, we look qualitatively at the effect of the AC 
model compared to DGP on various tests of the cosmological expansion history. 
Several cosmological datasets are typically used to constrain the expansion history at 
various epochs: type la supernovae , large scale structure , and the microwave 
background p, p^ . 

Type la supernovae are relatively reliable standard candles, and provide a good 
constraint on the recent expansion of the universe via the redshift-luminosity relation 
based on the luminosity distance di: 



(1 + ^) 
HQ\/\VLk 



(3.7) 



where S {X) = {X, sin X, sinh X) for a flat, closed or open universe respectively. Since 
the Hubble parameter is higher in AC cosmologies (at fixed flm) and increases with 



increasing f3, (|3.7| ) shows that this results in a lessening of the luminosity distance 
and hence a lower magnitude. Figure ^ demonstrates this with a direct redshift- 
magnitude plot. 

A more conventional visualization of the effect of the AC model is given by 
plotting the preferred regions of {Qm, a^} parameter space at different values of 
flk, P- Figure H shows the projection on the {f^m, oP'} plane at fixed (3 and fixed Vt^ 
values respectively: is plotted against VLm as this more readily compares with the 
VLr^ parameter conventionally used in the analysis of DGP models. The left figure 
indicates how the preferred region of parameter space reacts to the (3 parameter (two 
values, (3 = 0,a/2 are shown), and the right figure how the parameter space reacts 
to Qk ioT general (3. In each case the figure shows that decreases and a increases 
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0.0 0.5 1.0 1.5 2.0 

Redshift 

Figure 2: Plot of magnitude vs. redsliift for (from top to bottom curve) ACDM (green), 
DGP (blue) and the AC model with (q = 0.702, /? = 0.702) (red), along with the supernova 
redshift data. 0,^, = 0.27. 

in response to increasing (3, although it is interesting to note that the projection at 
fixed flk is relatively insensitive to that value of Qk- 

It is clearly not difficult to reproduce the supernova redshift luminosity relation in 
isolation, particularly if the possibility of an open universe is included. However, the 
real tension for DGP (and even more so for the asymmetric model) is in combining 
the supernova constraints with the constraints from other cosmological data . 
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Figure 3: An illustration of the restriction due to the Supernova data on {Jim; o^} param- 
eter space projected onto fixed fifc and fixed (3 subspaces. On the left, varies freely, and 
two fixed values of f3 are shown: 13 = the solid contour, and [3 = a/2 the dashed contour 
(the separate lines indicate Vtk = for each (3 value). In the right figure, [3 now varies freely 
and the bands indicate three fixed values of O^: from lightest to darkest = 0.1, 0, —0.1 
respectively. The separate lines indicate the bounding values of /?, the lower /3 = 0, and 
the upper (3 = a, where the lines are dotted for 0^, = 0.1, dashed for fi^ = 0, and solid 
for = —0.1. The large values of Vtk chosen indicate that in this projection of parameter 
space, the model is relatively insensitive to Vt^. 



The CMB shift parameter or (essentially) the ratio between the angular 
diameter distance to and horizon size at decoupling is typically used to constrain 



dark energy models |28|, as it is relatively model independent: 



-R(z* 



Vim Hi 



C 



1 + z^)DAiz^, 



/a 



dz' 



(3.8) 



where = 1090.51 ± 0.95 is the redshift at decoupling (and c has been tem- 
porarily reintroduced for reference). The problem with lowering Qm now becomes 
more apparent. While we can ensure that the comoving distance is maintained by 
dropping Q^n, the shift parameter is also lowered by this process. Indeed, fiat DGP 
requires fi^ ^ 0.35 to match the WMAP 5 year value R = 1.71 ± 0.02 g |28|. 



In order to compare the shift parameter constraint on the AC model to the 
situation with the DGP model we allow for open, fiat and closed cosmologies, 
and test the parameter space compatible with the given shift parameter. Figure 
^ shows allowed regions of {flm, ct"^} parameter space for three different (3 values 
ranging from the DGP to the Cardassian limit. These show that as (3 is increased, 
preferred values of Qm become higher, and plotting the allowed regions also shows 
how Qrn increases in response to increasing j3. Alternatively, we can take a different 
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projection by fixing and plot the allowed regions of {Qrn,C(^}, as indicated in 
the right hand figure of ^ In this plot, the limiting values of (3 are shown as lines, 
and increasing P corresponds to moving roughly diagonally upwards across the plot. 
Once again, this indicates that the preferred value of flm increases as /3 is increased. 



0.5 F 





0.1 0.2 0.3 0.4 0.5 0.6 

lOim 



0.1 0.2 0.3 0.4 0.5 0.6 



Figure 4: A depiction of the region of asymmetric cosmology parameter space consistent 
with the shift parameter. On the left, 0,^ varies freely, and three fixed values of (3 are 
shown: /5 = 0, or the DGP limit, is the lowest (blue) band, the green band an intermediate 
value of /?, and the grey band the maximal value of /?. The solid line indicates = 0. 
On the right, /? now varies freely and the bands indicate two fixed values of 0,k- the light 
band being a flat universe, and the dark band ilk = 0.03. The separate lines indicate the 
bounding values of /?, the lower /3 = 0, and the upper (3 = a, where the dotted lines are 
the limits for 0^ = 0, and the solid lines for $7^ = 0.03. 



We can now see how even just these two constraints on parameter space are 
problematic by combining them since increasing f5 tends to prefer a decreased Qm 
to fit the supernova data, yet an increased flm to fit the CMB shift parameter. In 
figure ^ we combine the plots, and include for reference an indication of the constraint 



coming from the baryon acoustic oscillation peak detected by the SDSS survey [§7 
This is usually represented as a dimensionless constant 







czi E{ziY/^ 



1 



dz' 



2/3 



= 0.469±0.017 
(3.9) 

where zi = 0.35, and Dy is the geometric average dilation scale |37|. There is some 
debate as to whether this measure should be used for models which do not behave 
as a constant equation of state dark energy |^, however, we include this band of 
parameter space as it seems likely that it is fairly indicative. 
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Figure 5: A look at the combined effect on the parameter space of asymmetric cosmology 
fixing Qfc (0 in the left figure and 0.03 in the right) and allowing (3 to vary between its two 
limits indicated by the lines. The dark (blue) band indicates the shift parameter preferred 
range of Vtm and a^, the lighter (grey) band from the supernova data, and the lightest 
(pink) band that from the BAO constraint. 

Figure ^ shows explicitly how increasing the (3 parameter steadily makes it more 
difficult to fit the various cosmological expansion data sets, as the regions of pa- 
rameter space consistent with near expansion history steadily diverge from those 
consistent with higher redshift data. We can see why this is the case by referring to 
the effective dark energy description. In spite of the enlarged parameter space, the 
asymmetric cosmology turns out to be under more cosmological tension than DGP. 
This is because for a given matter content, ( p.l|) shows that the Hubble parameter 
increases with redshift more rapidly than in DGP (which itself is more rapid than 
ACDM) as the j3 parameter increases. This means that for a given comovmg 
distance out to a particular redshift is lower in asymmetric gravity than DGP, which 
is correspondingly lower than ACDM. 



4. Discussion and Model Extensions 



There are other parameters one could include in both DGP and asymmetric cos- 
mologies. The general bulk spacetime of a cosmological braneworld includes a bulk 



black hole [33, 40], and while the effect of this black hole has been considered for 
Randall-Sundrum cosmologies (where it gives rise to a dark radiation term) it has 



not generally been included in DGP cosmologies (though see ||20| for a discussion of 
the problems it gives rise to for DGP). Adding in this general mass term as in ( ^.13| ) 
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alters (|]23D to 



^-''4{"'-'i)-'^<{M^,~k,) (4.1) 

where is the mass of the bulk black holes in the adS interior to the left of the 
brane, and /i2 an (effective) black hole mass of the exterior adS bulk on the RHS of 
the brane. 

Clearly, having a bulk black hole on the LHS (/^i > 0) simply adds in a 'dark 
radiation' term in the effective cosmological energy density in an analogous fashion 
to Randall- Sundrum cosmology [^. However, the effect of a black hole term on the 
RHS of the bulk is more interesting. Since the part of the bulk being excised on the 
RHS is the interior, we can have either sign for ^2 (see [^| for potential consistency 



problems with ^2 < 0), further, a positive mass bulk black hole actually leads to 
a negative contribution to the brane energy density. Setting /ii = 0, and writing 
= /^2/(-f^o'^o) fi^*^ that the effective Friedman equation is only subtly altered 
in the additional braneworld term: 

E{zf = a(l + zf^^^""^^ + 2a{a -(3) + 2a^{a - py + Q^{1 + 2)3(i+«'.) - Q^{l + z)^. 

(4.2) 

A negative fl^ (i.e. a negative black hole mass) simply adds to the value of E"^, and 
therefore will not assist the model in conforming to the expansion data. However, a 
positive bulk black hole mass contributes negatively, and therefore reduces the value 
of E"^. However, in order to prevent pressure singularities on the brane, we must 
ensure that < Qr, and thus the best that can be achieved by this term is a 
cancellation of the radiation density of the universe in the term under the square 
root, though not in the leading Einstein term. While this could lead to interesting 
effects in the early universe, these will be sub-leading and in any case it does not 
significantly help with fitting the late time expansion of the universe. 

To sum up: We have examined the asymmetric branes model p2|, ^ , a braneworld 
theory of modified gravity, with a view to exploring how well it can explain the late- 
time acceleration of the universe. The effective cosmological expansion above a Hub- 
ble distance of order 1mm is a one-parameter generalization of the DGP model, the 
effect of the extra parameter being to retard the expansion of the Universe relative 
to DGP. As such, it turns out that the asymmetric model has more problems fitting 
the cosmological expansion data than DGP. In addition, recent work on ghosts in the 



stealth model suggests that the AC model may well not be ghost-free around the 
accelerating vacuum, thus our overwhelming conclusion is unfortunately that pure 
AC models are not viable cosmological models for late time acceleration. Nonethe- 
less, it is important to check the behaviour of all possible concrete modified gravity 
models available to either identify or rule out alternatives to ACDM. 



- 14 - 



Acknowledgements 



We would like to thank Christos Charmousis, Anwar Gaungoo, and Antonio Padilla 
for useful discussions. EOC is supported by the European Commission's Frame- 
work Programme 6, through the Marie Curie Early Stage Training project MEST- 
CT-2005-021074, and AP is grateful to the University of Nottingham for financial 
support. 



A. Perturbation Theory and the Planck Mass 

Although the asymmetric model naturally lends itself to a boundary description 
of the equations of motion, in deriving the Planck mass it is useful to consider 
perturbation theory around the "domain wall" description, i.e. in which we take the 
full range of the coordinate y, and represent the brane as a physical delta function 
source: 



1 



Mf { Rab - -Rgab 



Perturbing these equations around the background (|2.9| ) yields: 



(A.l) 



T 



(A.3) 
(A.4) 



Following the procedure of Garriga and Tanaka |^2|], for constructing the Green's 
function, we see that the bulk solution for the spin 2 mode is /i^i, = Um{y)xt^\ where 
X is a 4D massive spin 2 tensor, and Um is found by solving ([A.2|) : 



AiJ2 



+ BiN2 



hi 



(A.5) 



where 



C = a-\z) 



Applying finiteness of the perturbation as y 
continuity and ( [A.2|) at the brane imply: 



(A.6) 

00 implies B2 = 0. Meanwhile, 



m 



Ml 



+ B1N1 



MIA2J1 



(A.7) 
(A.8) 
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Finally, normalization of the eigenfunctions gives 



lAr+i^ii 



m 



Thus our coefficients are completely specified as: 



Ai = Ao — 
2ki 



m 



-A 



where 



42 
^2 



4A;i 
TT^m 



/2 



Tim 
'2ki 



k2 
m 
ko 



J2 ( - 1 A^i ( - 1 - j^Ji (-]N2 



Ji 
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m 



Ml 



m 



Ml 
Ml 



Ji 



^ Mf' 



k2 
m 
ko 



J2 



(A.9) 

(A.IO) 
(A.ll) 



2M| 



Ml 



m 



kn 



m 



kn 



m 



Ji{^ J2{^ {N,l- N2[- +J,[- J2{t- 



m 



m 



m 



-1 



(A.12) 



Note that, as with the GRS and DGP models, there is no localizable zero mode, 
and the spectrum is continuous starting from = 0. 4D gravity must therefore be 
obtained as an effective behaviour within a range of scales. We therefore examine 
the Newtonian potential in the brane of a unit mass particle on the brane which is 
given in terms of the eigenfunctions by 



V{r) 



Mf 



'2 ^0 



dm- 



Airr 



\A2\'J? 



(A.13) 



Then, writing e = M|/Mj^, and Vc = 1/eki, and redefining the integration variable 
as a; = mrc/2 we have to leading order in s 



V{r) 



2 AkiE^ 
Ml 4nr 



dx- 



^g-(2r/rc)a: 



Ji(2x/k2rc) 
xJ2(2x/k2rc) 



(A. 14) 



4^4 



In order to estimate this integral, note that for r <^ r^, this integral is dominated by 
X ~ J1/J2 = where the integrand has a value of 0{s~^) with a width of 

Thus 



V{r) oc 



ki 1 



Ml 47rr 

and the Planck mass can be read off as m^^ ~ Ml/ki. 



(A.15) 
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